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ABSTRACT 

We present observations from our continuing studies of the earhest stages of massive 
star cluster evolution. In this paper, radio observations from the Very Large Array at 
0.7 cm, 1.3 cm, 2 cm, 3.6 cm, and 6 cm are used to map the radio spectral energy 
distributions and model the physical properties of the ultra-young embedded super star 
clusters in Henize 2-10. The 0.7 cm flux densities indicate that the young embedded 
star clusters that are powering the radio detected "ultradense H II regions" (UDH lis) 
have masses greater than ~ 10^ Mq. We model the radio spectral energy distributions 
as both constant density H II regions and H II regions with power-law electron density 
gradients. These models suggest the UDH lis have radii ranging between ~ 2 — 4 pc and 
average electron densities of ~ 10^ — 10^ cm~^ (with peak electron densities reaching 
values of ~ 10^ - 10^ cm ). The pressures implied by these densities are P/kB ~ 
10^ — lO^'' cm~^ K, several orders of magnitude higher than typical pressures in the 
Galactic ISM. The inferred H II masses in the UDH lis are ~ 2 — 8 X W^Mq; these 
values are < 5% of the embedded stellar masses, and anonamously low when compared 
to optically visible young clusters. We suggest that these low H II mass fractions may 
be a result of the extreme youth of these objects. 

Subject headings: galaxies: individual(Henize 2-10) — galaxies: star clusters — galaxies: 
star burst 
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1. INTRODUCTION 

The study of "super star clusters" (SSCs) has been revolutionized by the availability of high 
spatial resolution optical observations, such as those obtained with the Hubble Space Telescope 
(HST). These clusters appear to be common in starburst and merging galaxy systems (Whitmore 
2002), and theory suggests that extreme molecular gas pressures are a prerequisite for their for- 
mation (e.g. Elmcgrccn & Efrcmov 1997; Elmegreen 2002). Many SSCs have properties that are 
consistent with being adolescent globular clusters (although the necessary conditions for their sur- 
vival remain unclear, Gallagher &: Grebel 2002). 

Compared to globular clusters, SSCs are extremely young objects. However, from a formation 
perspective, everything has already happened by the time these massive clusters have emerged into 
optical light (Johnson 2002). Over the last few years, the study of SSCs has undergone a new 
revolution with the discovery of ultra-young SSCs that are still deeply embedded in their birth 
material. The pre-natal clusters have typically been identified in the radio (cm) regime as compact 
optically thick free-free sources having turnover frequencies > 5 Ghz. A sample of these objects 
have now been found in a number of galaxies (e.g. Kobulnicky &: Johnson 1999; Turner, Beck, &: 
Ho 2000; Tarchi et al. 2000; Neff & Ulvestad 2000; Johnson et al. 2001; Beck et al. 2002; Plante 
Sz Sauvage 2002; Johnson, Indebetouw, &; Pisano 2003). On a vastly smaller scale, objects with 
similar spectral morphologies are associated with extremely young massive stars in the Milky Way; 
these objects are known as "ultracompact H II regions" (UCH lis) (e.g. Wood &; Churchwcll 1989). 
UCH lis result from newly formed massive stars that are still embedded in their birth material and 
ionize compact (r << 1 pc) and dense (rig > 10"* cm^^) H II regions within their natal cocoons. The 
apparent similarity between UCH lis and the scaled-up extragalactic radio sources led Kobulnicky 
Sz Johnson (1999) to dub these extragalactic objects "ultradense H II regions" (UDH lis). 

Given the limited observations of UDH lis that are currently available, attempts to constrain 
their physical properties have been fairly crude. However, observations have provided estimates for 
the physical properties of these objects that are truly remarkable: their stellar masses can exceed 
~ 10^ Mq, the radii of their H II regions are typically a few parsecs, their electron densities suggest 
pressures of P/kB > 10^ K, and they appear to have ages of less than ~ 1 Myr (Kobulnicky 
& Johnson 1999; Vacca, Johnson, & Conti 2002; Johnson, Indebetouw, & Pisano 2003). Based 
on these estimates and comparisons to UCH lis in the Milky Way, a physical scenario for these 
ultra-young SSCs has been developed as that of compact (yet massive) clusters of stars that ionize 
extremely dense H II regions, which are in turn enveloped by warm dust cocoons. 

The dense H II regions can be observed in the radio regime, and the radio spectral energy 
distributions have turnover frequencies that increase with increasing electron density. The dust 
cocoons can be observed in the infrared to sub-millimeter and may have temperature profiles similar 
to individual massive protostars in the Milky Way (Vacca, Johnson, & Conti 2002), although the 
current lack of observations in these wavelength regimes does not allow rigorous constraints to 
be placed on the cocoon properties. For example, it is unclear whether the constituent stars are 
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surrounded by individual cocoons, or whether the cocoons have merged and the entire cluster is 
enveloped in a common cocoon. Presumably the relative morphology of the dust cocoon(s) and the 
stars depends on both the stellar density and evolutionary state of the cluster. Insufficient radio 
observations also make it difficult to model the properties of the embedded H II regions; while the 
nature of a compact radio source (non-thermal, thermal, or optically thick) can be determined from 
only a pair of high frequency data points, observations at a number of frequencies are required to 
accurately model spectral turnovers and constrain radii and densities. 

In order to distinguish between the different physical regions associated with pre-natal SSCs, 
throughout this paper we will use the following nomenclature: "cluster" or "SSC" refers only to the 
embedded stellar content, "UDH 11" refers to the dense H II region surrounding the stellar cluster, 
and "cocoon" refers to the dust cocoon surrounding the H II region. 

The starburst galaxy Henize 2-10 (He 2-10) is an excellent system in which to study the 
properties of ultra- young clusters in detail. At a distance of only 9 Mpc (Hq = 75 km s^^ Mpc~-^; 
Vacca k, Conti 1992), it is one of the nearest galaxies known to host multiple UDH lis (Kobulnicky 
&; Johnson 1999) as well as a large number of SSCs that have already emerged into optical and 
ultraviolet light (Johnson et al. 2000; Conti &; Vacca 1994). Mid-IR observations obtained by 
Vacca, Johnson, & Conti (2002) and Beck, Turner, & Gorjian (2001) confirm that the UDH lis 
are surrounded by warm dust cocoons. The dust cocoons surrounding the UDH lis in He 2-10 
are so luminous that Vacca, Johnson, &: Conti (2002) find they are responsible for at least 60% 
of the entire mid-IR fiux from this galaxy. This percentage of mid-IR flux due to clusters with 
ages < 1 Myr old is especially remarkable considering that He 2-10 has been undergoing an intense 
starburst for several Myr, and hosts nearly 80 optically visible SSCs (Johnson et al. 2000). 

The original discovery of the UDH lis in He 2-10 by Kobulnicky & Johnson (1999) utilized ob- 
servations at only two radio wavelengths in order to estimate their physical properties. In this paper 
we present radio observations of the UDH lis in He 2-10 at five wavelengths (6 cm, 3.6 cm, 2 cm, 

1.3 cm, and 0.7 cm) using relatively well-matched synthesized beams in order to more completely 
map their radio spectral energy distributions and better constrain their physical properties. 

2. OBSERVATIONS 

We obtained new Q-band (43 Ghz, 0.7 cm), K-band (22 Ghz, 1.3 cm), and U-band (15 Ghz, 
2 cm) observations of He 2-10 from February 2001 to January 2003 with the Very Large Array 
(VLA) ^. VLA archival data at U-band, X-band (8 Ghz, 3.6 cm), and C-band (5 Ghz, 6 cm) 
from May 1994 to January 1996 were also retrieved and re-analyzed. All of these data sets are 
summarized in Table 1. Based on the scatter in the VLA Flux Calibrator database, we estimate 
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the resulting flux density scale at each wavelength is uncertain by < 10%. 

Calibration was carried out in the Astronomical Image Processing System (AIPS) data re- 
duction package, including gain and phase calibration. The data sets at a given wavelength were 
combined, and the combined data sets were inverted and cleaned using the task imagr. While the 
uv-coverages at each wavelength are not identical, an attempt was made to obtain relatively well- 
matched synthesized beams by varying the weighting used in the imaging process via the "robust" 
parameter (robust= 5 invokes purely natural weighting, while a robust= —5 invokes purely uni- 
form weighting). The robust values used for each wavelength are listed in Table 2. The synthesized 
beams at each wavelength obtained with this method are relatively well-matched with the exception 
of the synthesized beam at 1.3 cm; the synthesized beam at 1.3 cm is approximately 1/3 the area 
of the beams at the other wavelengths. Therefore, the relative fluxes obtained at 1.3 cm should be 
regarded as lower limits. In order to facilitate comparison, the final images were all convolved to 
the same synthesized beam of 0."95 x 0."44 with a position angle of -3.2 degrees. Figure 1 shows 
the resulting 0.7 cm contours overlaid on the 3.6 cm grayscale. 

Flux densities of the sources were measured using two methods. In the first method, the task 
IMFIT in AIPS was used to fit the sources with two-dimensional gaussians; the uncertainty from 
these results was estimated by fitting each source using a range of allowed background estimates 
and Gaussian profiles. The second method utilized the VIEWER program in AIPS++ in order 
to place apertures around the sources at each wavelength; several combinations of apertures and 
annuli were used in order to estimate the uncertainty in this method. Both of these methods are 
fairly sensitive to the background level that is adopted, and the uncertainty in the resulting flux 
densities are dominated by this effect. The resulting flux densities and their uncertainties are listed 
in Table 3. With the exception of source #4, all of the 2 cm and 6 cm flux densities reported in 
this paper agree within uncertainty to the values presented by Kobulnicky & Johnson (1999); we 
attribute the single discrepant case to the different synthesized beams and background estimates 
used in each analysis. 



Table 1. VLA Observations of He 2-10 



A 


Antenna 


Date 


Obs. Time 


Flux 


Phase 


Phase Calib. 


(cm) 


Config. 


Observed 


(hours) 


Calib. 


Calib. 


(Jy) 


0.7 


C-array 


2003 Jan 05 


3.7 


3C286 


0836-202 


1.89 ± 0.06 


0.7 


C-array 


2002 Nov 21 


3.8 


3C286 


0836-202 


2.21 ± 0.04 


1.3 


BnA-array 


2001 Feb 04 


3.1 


3C286 


0836-202 


2.53 ± 0.05 


2.0 


BnA-array 


2001 Feb 04 


0.5 


3C286 


0836-202 


2.65 ±0.09 


2.0 


B- array 


1996 Jan 03 


2.5 


3C286 


0836-202 


2.98 ±0.01 


3.6 


B- array 


1996 Jan 03 


0.7 


3C286 


0836-202 


3.24 ± 0.01 


3.6 


A-array 


1995 Jun 30 


0.5 


3C48 


0836-202 


2.41 ± 0.06 


3.6 


BnA-array 


1994 May 14 


0.6 


3C48 


0836-202 


1.90 ± 0.01 


6.0 


A-array 


1995 Jun 30 


1.6 


3C48 


0836-202 


2.32 ±0.01 
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Fig. 1. — VLA 0.7 cm contours of He 2-10 overlaid on 3.6 cm grayscale. The contour levels are 
3, 5.0, 7, 8, 9, 10, 11 X a (0.07 mJy/beam). 
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3. PROPERTIES OF THE UDH lis 

3.1. Production Rate of Ionizing Photons 

In order to estimate the embedded stellar content of the UDH lis in He 2-10, their radio 
luminosities can be used to determine the production rate of Lyman continuum photons. Following 
Condon (1992), 

Q.,e > 6.3 X 10^^ s- (^) (-^) °" ( . (1) 

VlO^K/ VGhz/ V 1027 erg s-i Hz"^>' 

In the application of this equation (which assumes the emission is purely thermal and optically 
thin), it is advantageous to use measurements made at the highest radio frequency available for 
two reasons: (1) the higher the frequency, the less likely it is to contain a significant amount of non- 
thermal contaminating flux, and (2) the higher frequency emission suffers from less self-absorption 
and is therefore more likely to be optically thin. For these reasons, we use the 43 Ghz (0.7 cm) 
observations to determine QLyc for the UDH lis in He 2-10. An electron temperature must also 
be assTimed, and we adopt a "typical" H II region temperature of = 10^ K; the uncertainty in 
QLyc due to this assumption is < 20%. The resulting Qiyc values determined using this method 
are listed in Table 4, and they range from QLyc^7 X lO^^ - 26 x 10^^ s'^. To put these values 
in context, a "typical" 0-star (hereafter O*, the equivalent to type 07.5V; Vacca 1994; Vacca, 
Garmany, &; ShuU 1996) has an ionizing flux of QLyc = 1 x 10^^ s~^. Therefore, the QLyc values 
determined for the UDH lis in He 2-10 imply the equivalent of 700 - 2600 0*-stars in each cluster. 

There is an important caveat in regard to the QLyc values discussed above. A number of 
Galactic UCH lis are observed to be associated with diffuse extended emission. For example, in the 
Kim &; Koo (2001) sample of 16 UCH lis, all of the objects were associated with extended radio 
recombination line emission. Likewise, Kurtz et al. (1999) find that 12 out of 15 UCH lis in their 
sample arc associated with extended emission. These studies conclude that typically > 80% of the 
ionizing flux from the embedded stars in UCH lis is escaping to the outer envelope, possible due 
to a clumpy density structure within the parent molecular cloud. In the case of He 2-10, there also 
appears to be a somewhat diffuse thermal background in the immediate vicinity of the UDH lis. If 



Table 2. Imaging Parameters 



A 
(cm) 


Weighting 
(robust value) 


Synth. Beam 
("x") 


P.A. 

n 


RMS noise 
(mjy/beam) 


0.7 


1.0 


0.95 X 0.44 


-3 


0.07 


1.3 


5.0 


0.32 X 0.28 


-3 


0.05 


2.0 


1.2 


0.74 X 0.44 


-3 


0.05 


3.6 


0.8 


0.76 X 0.44 


-3 


0.03 


6.0 


2.0 


0.95 X 0.44 


-3 


0.05 
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Table 3. Radio Flux Densities of He 2-10 Sources 





Fo.7cm 


Fl.Scm 


F2cm 


Fs.ecm 


-^6 cm 


Source 


(mJy) 


(mJy) 


(mJy) 


(mJy) 


(mJy) 


1 


0.79 ± 0.13 


0.66 ±0.15 


0.91 ±0.12 


0.70 ±0.14 


0.61 ±0.18 


2 


1.08 ±0.14 


0.70 ±0.19 


0.97 ±0.14 


0.87 ±0.19 


0.66 ±0.17 


3 


0.23 ±0.05 


0.43 ±0.15 


0.68 ±0.13 


0.78 ±0.22 


0.66 ±0.18 


4 


2.91 ±0.57 


2.36 ± 0.48 


3.11 ± 0.38 


3.03 ±0.52 


2.26 ±0.36 


5 


0.89 ±0.21 


0.65 ±0.16 


0.83 ±0.15 


0.76 ±0.10 


0.58 ± 0.07 



Table 4. Estimated Properties of the UDH lis 









Constant Density Models 




Density Gradient Models 


# 


Q Lyc 


^'I stars 

(xW-'Mq) 


r 

(pc) 


Tie 

(xlO-^cm^^) 


Mhii 
(xIO^^Mq) 


7 


1-1/2 

(pc) 


"1/2 

(xlO'^cm^^) 


Mhii 
(xlO^M©) 


1 


7.1 


1.2 


-^•°-0.4 


7+3.2 


4.0 


-0.9 


1.7 


3.6 


5.5 


2 


9.7 


1.6 


l-°-0.4 


7 Q+3-7 
'••^-3.0 


4.4 


-1.3 


2.1 


2.4 


7.9 


4 


25.9 


4.3 


•J Q+1-2 
"^•^-0.9 




2.4 


-0.0 


3.1 


3.9 


2.4 


5 


7.9 


1.3 


1 7+0.4 
-^•'-0.2 


7 Q+1-6 
'••'-2.1 


3.7 


-1.1 


1.9 


2.9 


6.2 



^Qhyc values as determined from the 0.7 cm flux densities, assuming an HII region temperature of 10** K. 
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UDH lis are leaking ionizing flux in similar proportions to their Galactic counterparts, the inferred 
stellar content from the QLyc values determined above may be a significant underestimate. To 
estimate an upper limit on the possible leakage of ionizing flux from the UDH lis, the ionizing flux 
was measured from the entire region surrounding the UDH lis. The entire region has an ionizing 
flux of Qiyc ~ 8 X 10^^ s~^, approximately 40% higher than the sum of the ionizing fluxes from 
the individual UDH lis. This "extra" flux is only an upper limit on the ionizing flux that could 
be escaping from the UDH lis as it could also be due to contributions from other more evolved 
clusters in the region. 

The masses of the embedded stellar clusters can be estimated from their Lyman continuum 
fluxes using the Starburst99 models of (Leitherer et al. 1999). For a cluster less than 1 Myr old 
formed in an instantaneous burst with a Salpeter IMF, 100 Mq upper cutoff, 1 Mq lower cutoff 
(note that reducing the lower mass cutoff to 0.1 Mq increases the cluster mass by a factor of ~ 2.5), 
and solar metallicity, a lO^M© cluster has Qiyc ~ 6 x 10^^ s"-*^. Assuming that Qhyc scales directly 
with the cluster mass, we flnd that the UDH lis in He 2-10 are powered by stellar clusters with 
masses ~ 1.1 — 4.3 x 10^ Mq (Table 4). The mass derived for source 4 from these radio observations 
is in good agreement with, but a factor of 2 lower than the mass derived by Vacca, Johnson, & 
Conti (2002) using the bolometric luminosity of the cluster (after taking into account the 0.1 Mq 
lower mass cutoff used by Vacca et al.). While the small difference between the mass derived from 
the radio and bolometric luminosities is certainly contained within the uncertainties, the direction 
of the difference possibly suggests that a fraction of the ionizing photons are being absorbed by 
dust within the UDH II. Given the possible leakage of ionizing photons from the UDH lis and/or 
their absorption by dust, these mass estimates derived from the radio luminosities are lower limits, 
and the cluster masses could be as much as ~ lOx higher than these values. 



3.2. Comparison to Model H II Regions 

In order to gain additional physical insight about these H II regions, we can invoke simple 
models. Following the models of Johnson, Indebetouw, k. Pisano (2003), we model the compact 
radio sources as spherical H II regions with constant electron temperature of 10^ K under two 
different assumptions: (1) the electron density proflle is constant, and (2) the electron density 
profile is allowed to vary as a power-law of the form rig oc . For both sets of models, the flux 
density measurements at 43 Ghz, 15 Ghz, 8 Ghz, and 5 Ghz were used to fit the models (in a 
least-squared sense); the 22 Ghz data points were excluded from the models fits as they are only 
lower limits. 

The major difference between these models and those used by Kobulnicky & Johnson (1999) 
and Johnson et al. (2001) is the allowance for a power-law density profile within a given H II region. 
We also include limb-darkening, which has an effect in the optically-thin regime at high frequencies. 
However, the results we obtain with these new models are in excellent agreement with the results 
from the models used in the above papers; this provides additional confidence in the robustness of 
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the model results. However, the additional radio data points presented in this paper allow us to 
place more rigorous constraints on the model parameters than Kobulnicky & Johnson (1999). 

For the models assuming constant electron density, the data are best fit by H II regions with 
radii that range between r« 2 — 4 pc and electron densities of ng ~ 4 — 7 x 10'^ cm^^. These best 
fit models are shown with the dashed line in Figure 2. The modeled radii and densities were then 
used to determine total H II masses of the UDH lis of ~ 4 — 24 x 10^ Mq. Given that we have two 
free parameters (n-e and r) and four data points, the models are over-constrained and we determine 
the possible range in the model parameters by finding the minimum and maximum values for rig 
and r that can fit the data within uncertainty. All of these properties are listed in Table 4. 

The second set of models allowed for a density gradient of the form ng oc r'"'. The main affect 
of such a gradient on the modeled spectral energy distributions is to make the turnover less abrupt. 
For sources 1, 2, and 5 these models produce a marginally better fit than the constant density 
models; the best fit model for source 4 had 7 = 0, which is equivalent to the constant density 
model. These best fit models are shown with a solid line in Figure 2. For these models, the best 
fit density gradient (7), the half-mass radius (ri/2)) the electron density at the half-mass radius 
{ni/2), and the total inferred H II mass (Mhii) are listed in Table 4. For the sources best fit with 
a density gradient, the electron densities exceed values of f« 10^ cm~^ in the inner 0.5 pc, and the 
core densities can reach values close to ~ 10^ cm~^. 

The inferred H II mass for each UDH II is < 5% of the embedded stellar mass (~ 3% -5% 
for sources 1, 2, and 5, and ^ 1% for source 4), which is a strikingly low value when compared 
to optically visible young clusters. For example, the R136 cluster in the 30 Doradus nebula has a 
stellar mass of ~ 6 x 10^ Mq (Hunter et al. 1995), but the H II associated with the region has a 
mass of ^ 1 x 10^ (Peck et al. 1997). As the radio observations presented here do not suffer from 
extinction, it is not likely that the low H II mass fraction for these objects is due to intervening 
absorption. Instead, we suggest that the low H II mass fraction might be due to the extreme youth 
of the stellar clusters; much of the gas mass in the vicinity of the cluster may still be shielded from 
the ionizing radiation, and therefore still in molecular or neutral atomic form. 

The densities suggested by the models described above imply tremendously high pressures of 
P/ks ~ 10^ — lO^'' cm~^ K within the UDH lis. Pressures of this magnitude are also typical 
in Galactic UGH lis (e.g. Ghurchwell 1999, and references therein), but they are extremely high 
compared to typical pressures in the ISM oiP/kB ~ 10"^ — 10^ cm~^ K (Jenkins, Jura, &: Loewenstein 
1983). These high pressures may be one of the requirements for the formation of SSGs (e.g. 
Elmegreen & Efremov 1997; Elmegreen 2002). However, there are two caveats to bare in mind 
when considering these estimated pressures: (1) A cluster itself will contribute to the pressure of 
the surrounding medium. We can estimate the affect of the cluster by assuming the natal H I 
cloud had a gas temperature of Thi ^ lO^iC, and the photoionization of the cluster increases the 
temperature to Thii ~ lO^K; therefore, the cluster itself will cause the pressure to increase a factor 
of ~ 10^ over the initial value in the natal cloud. (2) We are not directly measuring the initial 
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Source #1 Source #2 




Source #4 Source #5 




10 20 30 40 50 10 20 30 40 

V (Ghz) V (Ghz) 



Fig. 2. — H II region model fits to the UDH lis in He 2-10 (source 3 is not included as its spectral 

energy distribution cannot be fit by a thermal model). The dashed lines are the best fits from the 
constant density models, and the solid lines are the best fits from the models allowing for a density 
gradient. Because the synthesized beam at 22 Ghz is significantly smaller than the beams at the 
other wavelengths, the flux density measured at 22 Ghz is a lower limit and was not used to fit the 
models. 
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pressure of the newborn H II region. Despite the youth of these objects, the H II regions have 
had some time to expand toward pressure equihbrium. Consequently, the pressures inferred at this 
stage in their evolution are lower than the initial values. If an H II region expands at a sound 
speed of roughly ~ 10 km/s, it will have expanded from ~ 1 pc to ~ 10 pc in 1 Myr, causing the 
pressure to decrease by roughly a factor of ~ 10^. These are fairly crude estimates, however even 
in the case that effect (1) is dominant, the pressures implied for the quiescent natal material are 
P/kB ~ 10^ - 10^ cm-=^ K. 

Source 3 cannot be fit by any purely thermal model (Figure 3). Although this source appears 
to have a slight turnover around 5 Ghz, the data at the remaining frequencies are fit by a fairly 
steep negative power-law of a ~ — 1 (where Si, oc v"), indicating a non-thermal origin. The nature 
of this source is discussed in Vacca et al. (in prep). 

4. SUMMARY 

Multi-frequency radio observations of the deeply embedded massive star forming regions in 
He 2-10 have allowed us to better constrain their physical properties. The ionizing fluxes calculated 
from the thermal radio continuum have values ranging from QLyc ~ 7 x 10^^ — 26 x 10^^ s~^, which 
imply the clusters have stellar masses > IO^Mq. These masses are fully consistent with the masses 
estimated for the optically visible SSCs by Johnson et al. (2000). We model the H II regions under 
the assumptions of both constant electron densities and a power-law electron density gradients. 
These models suggest that the dense H II regions have radii < 4 pc (and as small as 1.7 pc), 
and the electron densities may reach values as high as ~ 10^ cm~^ in the cores of these regions. 
The inferred H II masses for the UDH lis are < 5% of the embedded stellar masses, anomalously 
low compared to optically visible young clusters, and possibly due to the extreme youth of these 
objects. The densities derived from the models imply pressures of P/kB ~ 10^ - 10^° cm K 
within the UDH lis, typical of Galactic UCH lis, and provide us with observational confirmation 
of the extremely high pressures involved in the early stages of super star cluster evolution. 

We are grateful the NRAO staff for their assistance, and Claire Chandler in particular for her 
advice on the high frequency observations. We thank Bill Vacca, Miller Goss, and Jay Gallagher 
for many useful discussions on this subject and their comments on the manuscript. Remy Indebe- 
touw kindly provided assistance developing the physical models presented in this paper. We also 
thank the anonymous referee for providing stimulating comments that led to improvements in the 
manuscript. K.E.J, gratefully acknowledges support for this paper provided by NSF through an 
Astronomy and Astrophysics Postdoctoral Fellowship. 
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Source #3 




Fig. 3. — The spectral energy distribution of source 3. This spectral energy distribution cannot be 
fit by a thermal model. 
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